Power Laws in Solar Flares: Self-Organized Criticality or Turbulence? 
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We study the time evolution of Solar Flares activity by looking at the statistics of quiescent times 
tl between successive bursts. The analysis of 20 years of data reveals a power law distribution with 
exponent a ~ 2.4 which is an indication of complex dynamics with long correlation times. The 
observed scaling behavior is in contradiction with the Self-Organized Criticality models of Solar 
Flares which predict Poisson-like statistics. Chaotic models, including the destabilization of the 
laminar phases and subsequent restabilization due to nonlinear dynamics, are able to reproduce the 
power law for the quiescent times. In the case of the more realistic Shell Model of MHD turbulence 
we are able to reproduce all the observed distributions. 

PACS Number(s): 96.60.Rd; 47.52+J; 05.65+b 
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Solar flares are sudden, transient energy release above 
active regions of the sun Q. Energy is released in vari- 
ous form (thermal soft X-ray emission, accelerated par- 
ticles, hard X~ray (HXR) emission, and so on). Parker 
1^ conjectured that flares represent the dissipation at the 
many tangential discontinuities arising spontaneously in 
the bipolar fields of the active regions of the Sun as a con- 
sequence of random continuous motion of the footpoints 
of the field in the photospheric convection (2| . Probabil- 
ity distributions, calculated for various observed quan- 
tities, X, can be well represented by power laws of the 
form P{x) = Ax~°' . In particular from HXR emission, 
the distribution of peak flux yields a ~ 1.7, that of total 
energy associated with a single event yields a ~ 1.5 and 
finally the distribution of flare duration yields a ~ 2 . 

The conjecture by Parker and the power laws found 
in the distributions of real flares stimulated a new way 
of looking at impulsive events like flares. In fact Lu and 
Hamilton [Q pointed out that Self-Organized Criticality 
(SOC), introduced earlier could describe the main 
features of HXR flares, [^j7j, even if recent 2D MHD sim- 
ulations have been devoted to recover power laws in the 
energy dissipation What is usually called SOC is a 
mechanism of charging and discharging, apparently with- 
out tuning parameters, which reproduces self-similarity 
in critical phenomena. 

Lu and Hamilton, [Q, assume that the coronal mag- 
netic field evolves in a self-organized critical state in 
which an "event" can give rise to other similar "events" 
through an avalanche process (sandpile model). An ac- 
tive region on the sun is thus modeled through a magnetic 
field B on a uniform 3D lattice. In order to have a statis- 
tically stationary state, energy is injected into the system 
by adding a small magnetic field increment SB at a ran- 



dom site on the grid. When an avalanche takes place, the 
energy input is suspended until all sites become stable. 
In this sense the avalanches (flares) are fast phenomena, 
on a time scale much smaller than the injection mech- 
anism. The large popularity of these models settles on 
their capability to reproduce the power law behavior in 
the distribution functions of the total energy, the peak 
luminosity and the duration of avalanches. 

What we want to stress in this letter is the fact that 
a different kind of statistics can be studied on solar flare 
signals: the distribution of laminar or waiting times, 
i.e. the time intervals between two successive bursts. 
This distribution has been recently studied on solar flares 
HXR events |Jl^. Wheatland et al. have also em- 
phasized the fact that this kind of distribution is crucial 
from the point of view of the avalanche model. SOC mod- 
els indeed are expected to display an exponential waiting 
time distribution P{tl) = (r^)^^ exp("TL/(TL)), where 
{tl), the average laminar time, depends on the param- 
eters of the model. This behavior is related to the fact 
that the avalanche duration is much smaller than the 
charging time (the time between two successive throws 
of magnetic field in random position) and charging place 
is independent on the avalanche position. Then one ex- 
pects no correlation between successive bursts and thus 
a trivial statistics for the laminar times. This is clearly 
observed in a simulation of the SOC automaton that we 
have done (see the inset of figure 0). On the contrary 
all authors |^,|l^ found a more or less well defined power 
law distribution. In particular Wheatland et al. [|0|, by 
performing a careful statistical analysis of waiting time 
distribution on 8 years of solar flares HXR bursts ob- 
served by the ICE/ISEE 3 spacecraft, have shown that 
the distribution in no way can be attributed to a nonsta- 
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tionary Poisson process. 

We have done the same statistics on laminar phases 
using twenty years of data from National Geophysical 
Data Center of USA. In this database starting and end- 
ing times as well as peak times of HXR bursts associated 
to flares from 1976 up to 1996, measured at the Earth by 
satellites in the 0.1 to 0.8 nm band, are stored. We calcu- 
lated the laminar times as the time differences between 
two successive maxima of the flares intensity recorded 
during periods of activity of the same instrument. Two 
different kind of analysis have been performed: in the 
first one we have built up a dataset of about 1100 sam- 
ples (hence on dataset A) by calculating only the differ- 
ences between the time of occurrence of flares within the 
same active region, as identified through the Ha flares 
occurrence. To build up the second dataset (dataset B), 
we have considered the sun as a unique physical system, 
and we have calculated the time differences between two 
successive maxima of flare intensity regardless of the po- 
sition of the flare on the sun surface. In this way we 
get a dataset of about 32, 000 samples. The analysis of 
these data shows that, in both cases, laminar times dis- 
play a clear power-law distribution P{tl) = AtJ^°' with 
a = 2.38 ± 0.03 in the range 6 hr < < 67 hr (re- 
duced = 2.2) for dataset B and a = 2.4 ±0.1 (reduced 
= 1.1) for dataset A (figure Q). The exact value of 
the exponent can be affected by the finite lenght of the 
observation times, which underestimates the occurrence 
of long waiting times. However these results, as well as 
those obtained by previous authors ]s|JlO[|, allow us to 
consider the power law distribution of waiting times as 
firmly established as the power laws observed for total 
energy, peak luminosity and time duration and force us 
to investigate whether models different from SOC can 
account for all these distributions. 

One point of the SOC philosophy is that the system has 
to be at the edge of chaos in order to display power laws 
. Indeed in a chaotic system there exists a character- 
istic time tc given in terms of the leading Lyapunov ex- 
ponent as tc — 1/A. This seems to be incompatible with 
the existence of self-similar scaling laws (i.e. no charac- 
teristic times) a part the limit case A = 0. The above 
argument in not very strong. As a matter of fact we re- 
call that it is possible to have chaotic systems showing 
scaling behavior in presence of many time scales or in 
presence of strong fluctuations of the local Lyapunov ex- 
ponent [|l3| . Thus it is worth investigating which kind of 
chaotic systems can give rise to self-similar scaling laws. 

It is worth noting that the occurrence of a power law in 
the distribution of the laminar times is the analogous for 
the Solar flares of the Omori's law for the earthquakes 
[0 and represents a clear indication of the existence, 
in the flare dynamics of strong correlations between suc- 
cessive bursts, at variance with the SOC model. The 
unique possible origin of the correlations arises from non 
trivial evolution equations of the phenomenon. Because 
solar flares are governed by MHD equations, it is rather 
natural to investigate their statistics in the context of 



MHD turbulence. Turbulence is a common phenomenon 
in fluids, where chaotic dynamics and power law statis- 
tics coexist. Moreover fluid turbulence displays time in- 
termittency, in that dissipative events are not uniformly 
but burstly distributed in time. 

Before introducing the turbulence model, let us stress 
the fact that the intermittent behavior can be observed 
also in simple dynamical models. Let us consider the one 
dimensional random map xt+i = rtXt{\ — xt), where at 
each step t the random variable rt is extracted according 
to a given distribution, e.g. = 4 with probability p and 
rt — 1/2 with probability 1 — p. If p < 1/3 a; = is an 
attracting fixed point. If p — 1/S + 6p one has a rather 
interesting behavior called on-off intermittency |T5| : Xt 
remains close to zero for a certain time (laminar time) 
then there is a short interval of strong activity (burst) 
after which a new quiescent phase takes place. Also in 
this simple model, if Sp is not too large, one observes 
power law not only for the energy distribution (here de- 
fined as the integral of Xt over the burst), but also for 
the laminar times r^. The exponent of power law for tl 
turns out to be a = 1.5 [l^ . In spite of its simplicity, the 
random map model contains some basic ingredients of the 
relevant features of the time intermittency in dynamical 
systems, i.e. the destabilization of the laminar phase by 
linear instability and the subsequent restabilization due 
to the nonlinear dynamics. 

Dynamical systems more directly related to fiuid tur- 
bulence are the so called shell models [|l^. Shell models 
represent a zero-order approximation of fluid equations 
(Navier-Stokes or MHD equations) in which one consider 
a single (complex) scalar variable u„ (and bn) as rep- 
resentative of the velocity (magnetic) fluctuation associ- 
ated to a wavenumber A;„ — fco2" {n = 1, N). The fact 
that the wavenumbers kn are exponentially spaced allows 
to reach very large Reynolds numbers with a moderate 
number of degrees of freedom and then to investigate 
regimes of 3-D MHD turbulence which are not accessible 
by direct numerical simulation. The philosophy underly- 
ing the shell model approach to turbulence is that even 
with a relatively small dynamical system it is possible to 
reproduce some statistical features of the turbulent cas- 
cade. In particular shell models mimic at the best the 
time intermittency of real turbulent fluid flows. 

The evolution equations for the dynamical variables Un 
and bn are built up by ignoring any detail of the spatial 
structure and boundary conditions. Only the interac- 
tions between nearest and next nearest neighbor shells 
are retained in the form of quadratic nonlinearities. The 
coupling coefficients of nonlinear terms are determined by 
imposing the inviscid conservation of the MHD quadratic 
invariants jl^Jl^. The particular shell model we used in 
our simulation reads ]T9[ | 

= -lyklun + fn+ ikn {(m„+iU„+2 - bn+lbn+2) " 
- {Un-lUn+1 - &„-l6,i+l) - - {Un^2Un-l ^ fen-2&ri-l)}* 
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= -rjk^hn + ifc„(l/6) {(u„+ifo„+2 - &«+iu«+2) + 

(m„_i6„+1 - 6„_iU„+i) + {Un-2bn-l - 6ri-2W«-l)}* 

where and 77 are respectively the viscosity and the re- 
sistivity and /„ is an external forcing term acting only 
on velocity fluctuations. 

Shell models are good models of turbulent cascade in 
the sense that they display, in the limit of fully developed 
turbulence i^, 77 — > 0, scaling laws for the structure func- 
tions Sp{n) = (|a;„|P) fc„^'' where a;„ is either u„ or 6„. 
In the hydrodynamic limit (&„ = 0) the set of scaling ex- 
ponents Cp ai'e found to be very close to those obtained by 
experiments One indeed observes a clear deviation 
from the Kolmogorov scaling {(^p — p/3) as a consequence 
of the intermittent dynamics of the system. 

Another observable whose statistics is well reproduced 
in shell models is the energy dissipation e{t) defined as 

= '^E^^i^nl'^nP + vT,n=i'^lK\'^ which displays 
the characteristic intermittency of fully developed turbu- 
lence (see figure ||). 

We have performed on e{t) the same statistical anal- 
ysis done for the solar flare signal. We deflne a burst 
of dissipation (corresponding to a flare) by the condition 
e(^) > Cc- This definition allows us to calculate the distri- 
bution functions for the peak values of the bursts, their 
total energy (defined as the integral of the signal above 
Ec) and the duration of the bursts (defined as the time 
during which the dissipation is above ec). We have cho- 
sen the threshold as Cc = ( e(i) ) -I- 2a, where the average 
and the standard deviation have been calculated on the 
time intervals in between the bursts, through an iterative 
process in order to take into account only the background 
contribution. The results of the analysis are shown in fig- 
ure ^ Also in this case we observe clear power law dis- 
tribution functions with exponents a ~ 2.05 for the peak 
distribution, a ~ 1.8 for the total energy distribution and 
a ~ 2.2 for the burst durations. The exponents are close 
to those obtained in analyzing solar flares data, but we 
do not think that the agreement is particularly signifi- 
cant as the model exponents depend on the value chosen 
for the threshold. The relevant point is that, at variance 
with SOC models, MHD shell models display a power law 
statistics also for the laminar times, as shown in figure ^. 
The scaling exponent turns out to be a ~ 2.70, close to 
the one obtained from the experimental data. 

The different behavior of SOC models and turbulent 
MHD shell models is related to the conceptually differ- 
ent mechanisms underlying the SOC phenomenon and 
the phenomenon of intermittency in fully developed tur- 
bulence. SOC models represent self-similar phenomena, 
while the intermittent behavior of turbulence is related 
to its chaotic nature Let us stress again that the 

actual value of the numerical scaling exponent is not im- 
portant as it could depend on the details of the model. 
The statistics of the laminar times between two bursts is 
due to global properties of the system and it is thus more 
relevant than the properties of the single burst. In this 



sense, a good model for the flare bursts should be able to 
reproduce the quiescent time distributions. SOC models 
predict a Poissonian statistics for the scaling behavior of 
the quiescent times in the Solar flares activity. On the 
contrary chaotic models are able to reproduce the power 
law for the quiescent time since they include the correct 
mechanism for destabilization of the laminar phases and 
subsequent nonlinear restabilization. 
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FIG. 1. Probability distribution of the laminar time P(ri) 
between two X-ray flares for dataset A (dashed line) and 
dataset B (full line). The straight lines axe the respective 
power law fits. In the inset we show, in lin-Iog scale, the 
distribution for dataset B (full line) and the distribution ob- 
tained through the SOC model (dashed line) which displays a 
cleax exponential law. The variables shown in the inset have 
been normalized to the respective root-mean-square values. 
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FIG. 2. Time series of energy dissipation e{t) for the shell 
model. The parameters used in the simulation arc A'" — 19, 
u = rj = 10~^, fco = 1. The external forcing term /„ is a 
stochastic variable acting only on the first two shells of the ve- 
locity fluctuations. It is calculated according to the Langevin 
equation dfn/dt = —fn/ro + n, where To is the characteristic 
time of the largest shells (ro ~ 10 in our units) and fi is a, 
Gaussian white-noise with a = 0.1. 
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FIG. 3. Total energy distribution P(e), energy peak dis- 
tribution P{p) and bursts duration distribution P{tb) for the 
shell model. The variables have been normalized to the re- 
spective root-mean-square values. The straight lines are the 
fits with power laws. The values of P(e) and P{tb) are offset 
by a factor 100 and 10^^ respectively. 
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FIG. 4. The distribution of laminar times P{tl) for 
the Shell model, normalized to the root-mean-square. The 
straight line is the fit with a power law. 
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